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Abstract 

We report a search for production of long-lived charged massive 
particles in a data sample of 90 pb“^ of y/s = l.STeVpp collisions 
recorded by the Collider Detector at Fermilab (CDF). The search uses 
the muon-like penetration and anomalously high ionization energy loss 
signature expected for such a particle to discriminate it from back¬ 
grounds. The data is found to agree with background expectations, 
and cross section limits of 0(l)pb are derived using two reference 
models, a stable quark and a stable scalar lepton. 

PACS Numbers: 13.85.Rm 12.60.Jv 14.80.Ly 

Many models for new physics introduce new particles which can be long- 
lived either due to a new conserved quantum number (e.g., R-parity in su¬ 
persymmetry) or because the decays are suppressed by kinematics or cou¬ 
plings 1^0]. If they are electrically charged, these particles can be detected 
directly. The possibility of new charged particles which are absolutely sta¬ 
ble is constrained by cosmological considerations and by searches for exotic 
particles in stable matter [Q. However, particles which are not absolutely 
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stable but are long-lived on an experimental scale (100 ns) are constrained 
only by direct searches. The most stringent limits are set by a previous 
search at Fermilab’s Tevatron collider [Q and by searches at CERN’s LEP2 
collider 0 probing masses up to about 90 GeV/c^. In this letter, we present 
the results of a new search for production of long-lived charged massive par¬ 
ticles (CHAMPs) using a data sample of 90 pb“^ of y/s = l.STeVpp colli¬ 
sions recorded by the Collider Detector at Fermilab (CDF) during 1994-5. 
We search for particles with anomalously high ionization energy loss rate, 
dE/dx, which would be produced by a slow massive charged particle. 

The search can be applied to several models which fall naturally into 
two distinct categories; weakly produced particles (e.g., new leptons), and 
strongly produced particles (e.g., new quarks). The lower production cross 
section of weakly produced particles only yields sufficient events for masses 
.5100 GeV/where the background is high, while the higher cross section 
of strongly produced particles allows sensitivity at higher mass where the 
background is low. The search is made as model independent as possible, 
but to quantify the results we use a long-lived fourth generation quark as a 
reference model for a strong production search and Drell-Yan production of 
a long-lived slepton from gauge-mediated supersymmetry breaking (GMSB) 
scenarios for a weak production search. 

The CDF detector, described in detail in Ref. p, measures the trajec¬ 
tories (tracks) and transverse momenta [0, pt, of charged particles in the 
pseudorapidity region {pi < 1.1 with the central tracking chamber (CTC) 
and silicon vertex detector (SVX), which are immersed in a 1.4 T solenoidal 
magnetic held. Up to 54 time-over-threshold measurements made by the 
CTC for each track determine the dE/dx with an average resolution of 13%. 
The charge deposited in each of the four layers of the SVX provides a sec¬ 
ond measure of the dE/dx with an average resolution of 18% [Q. Control 
samples with well identihed particle types are used to calibrate the dE/dx 
measurements at different velocities; electrons and muons from W boson de¬ 
cay at high velocity (/^y > 100), muons from 3/%/ decay and pions Ks decay 
at intermediate velocity, and protons and deuterons from secondary interac¬ 
tions in the beampipe at low velocity (/^y < 1). Fig. shows the comparison 
of these measurements to the predictions. Electromagnetic and hadronic 
calorimeters, located outside the superconducting solenoid, measure energy 
in segmented p—cj) towers and identify electron candidates. Drift chambers for 
muon identihcation are situated outside the > 5.3 interaction lengths (Aint) 
thick calorimeters and behind an additional > 3.5Aint thick steel absorber. 
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Three different trigger data sets are used for this search. A muon trigger 
selects events with hits in the muon chambers which match a track with 
Pt > 12 GeV/c in the CTC within 5°. A massive particle can penetrate 
the calorimeters and pass the muon trigger even if it is strongly interacting 
because the energy lost in hadronic interactions with the relatively light 
nucleons is too small to initiate a shower. The CHAMP mass is > 100 
times the nucleon mass so the energy available in the center-of-mass frame 
falls below the threshold for single pion production [^]. Only triggers in 
the region \t]\ < 0.6 are used because at larger \ri\ timing requirements that 
assume (3 = 1 are used to reduce backgrounds from beam losses. 

The second trigger selects events with missing transverse energy (.^r) 
> 35 GeV, which can arise since the CHAMPs will penetrate the calorimeter 
without fully depositing their energy p[. This trigger also provides accep¬ 
tance for events containing neutrinos, as is possible in GMSB models where 
CHAMPs are produced along with neutrinos from the cascade decay of a 
heavier particle. An electron trigger, which selects events containing elec¬ 
tron candidates within the range |? 7 | < 1.1 and with Ex > 18 GeV, provides 
additional acceptance for these cascade decays, as does the muon trigger, 
since charged leptons may be produced in these decays as well. 

The search selects charged particle tracks with \p\ > 35 GeV/c and 
\t]\ < 1 which have sufficient hits in the GTG and SVX to reduce back¬ 
grounds from misreconstructed tracks. The 35 GeV/c momentum cut is 
chosen because for lower momentum a GHAMP in the mass range of interest 
(M > 100 GeV/c^) would be moving too slowly to be efficiently recon¬ 
structed. The SVX and GTG dE/dx measurements are each required to be 
larger than the values expected for a particle with (d'y = 0.85. In the region 
/dy < 0.85, dE/dx (x 1//?^ to a good approximation. That allows calculation 
of a measured mass, M^E/dx, from the dE/dx and the momentum. The mass 
resolution is measured to be 20% using a calibration sample of protons and 
deuterons. The search is performed for different assumed mass, M, between 
100 and 270 GeV/c^ with 10 GeV/c^ steps. At each step, M^E/dx is required 
to be > 0.6 X M, a 2a cut. When combined with the dE/dx cut, this pro¬ 
vides additional background rejection at lower momentum. To be considered 
in the weak production search, tracks must additionally pass an isolation cut 
requiring less than 4 GeV of calorimeter energy or total track pt within a 
cone of -|- |A0P = 0.4 around the track. 

Backgrounds arise from tracks for which the dE/dx measurement fluctu- 
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ated high or included extra ionization from an unreconstructed overlapping 
particle. To determine the background, we use a control sample which is iden¬ 
tical to the search sample but at lower momentum (20 < |p| < 35 GeV/c) 
where signal would not contribute. The fake rate, dehned as the fraction 
of tracks in the control sample with dE/dx measurements high enough to 
correspond to (3'^ < 0.85, is measured to be 0(10“^) for all the different trig¬ 
ger datasets described above. The momentum dependence of the fake rate 
within the control sample matches expectations which allows us to extrapo¬ 
late the fake rate to the high momentum signal region. The probability of 
a high fluctuation in the dE/dx distribution obtained from this fake rate is 
used to scale the number of candidate tracks, which pass all selections except 
the dE/dx requirement, to obtain background predictions of 12 ± 2 tracks in 
the muon trigger dataset and 63 ± 9 in the.^T trigger dataset. The expected 
mass distribution for fake tracks in the signal region is shown in Fig. ^ It is 
obtained by folding the momenta of the tracks into the dE/dx distribution 
from the control sample with the assumption that large values of dE/dx are 
due to high mass particles. In the data, 12 and 45 tracks pass all cuts for 
the muon and^T trigger datasets respectively. Their mass distribution, also 
in Fig. 1^, shows no signihcant excess over the predicted background. 

For the weak production search, the isolation cut reduces the background 
to 0.85 ± 0.25, 4.0 ± 2.8, and 0.7 ± 0.5 tracks in the muon, ]/]rp, and electron 
trigger datasets respectively. In the data, 0, 1, and 0 tracks are observed in 
these samples. 

The signal efficiencies are determined using Monte Carlo simulation pro¬ 
grams and control data samples. The muon trigger efficiency is 80.5 ± 3.0%, 
and the track selection efficiency is 51.3 ± 2.5%, dominated by acceptance 
in the SVX. The tracking efficiency decreases at low velocity, /dy < 0.4, due 
to drift time limits in the CTC track hnding algorithms. This is measured 
with a sample of deuterons which are produced from secondary interactions 
in the beampipe. The efficiencies of the cuts on the kinematic variables \ri\, 
\p\, (3'y (dE/dx) and mass are model dependent. To set generally applica¬ 
ble limits, we determine these efficiencies using easily quantihable reference 
models. For the strong production case we use a long-lived fourth genera¬ 
tion quark calculated with the Pythia Monte Carlo program [^. The total 
efficiency increases from 1.5 to 2.9% over the mass range 100 — 270 GeV/c^ 
for a charge |e quark (U) and from 0.8 to 1.6% for a charge —|e quark (D). 
The charge asymmetry of the efficiency arises from the light quark (u, d, 
s) contributions to the fragmentation; Us and Ud mesons are charged while 


only the Du meson is charged. Furthermore, although a massive quark would 
efficiently penetrate the calorimeters, the hadron containing it can undergo 
charge exchange from interactions in the calorimeter which replace the light 
quark in it, and a |e quark is more likely to remain in a charged hadron and 
be detected by the muon chambers. The efficiency for this depends on the s 
quark suppression which is taken to be 30% The uncertainty from this 


effect, estimated by taking half of the efficiency difference obtained if every 
hadron is assumed to interact, is 20% for q = \e and 13% for q = |e. Other 


systematic uncertainties are 4% for trigger efficiency, 5% for track selection, 
4% for luminosity, and 7% from the choice of CTEQ3M as the parton 
distribution function. The total systematic uncertainties on efficiency are 
23% and 17% for q = \e and q = respectively. 

Figure ^ shows the cross-section limits we derive as a function of mass. 
From comparison with the expected cross-section, we derive mass limits at 
95% conhdence level of M > 190 GeV/c^ for q = \e and M > 220 GeV/c^ 
for q = |e. The charge exchange effects described above could be different 
for other models. To ease comparison with other models, we include in 
Fig. ^ a limit calculated without these effects. These limits are based on data 
collected with the muon trigger. The.^r trigger dataset is also searched since 
it could provide sensitivity to signal, but the trigger efficiency depends 
critically on the calorimeter’s response to a GHAMP which is very uncertain. 
This makes any cross-section calculations unreliable, so the^T trigger dataset 
is not included in the limit calculation for the strong production search. 

For the weak production search, the muon trigger and track quality cut 
efficiencies are similar to the strongly interacting case. The efficiencies of the 
model dependent kinematic cuts are estimated using as a reference model 
the Drell-Yan pair-production of stable sleptons calculated with the SPythia 

r4|. The total efficiency varies from 2.5% to 4.5% 
- 120 GeV/c^. The systematic uncertainties on these 


Monte Garlo program 
over the mass range 80 
efficiencies are similar to the strongly interacting case, without the charge 
exchange uncertainty. The cross-section limits obtained for direct slepton 
production range from 1.3 pb at M = 80 GeV/c^ to 0.75 pb at 120 GeV/c^ 

. The expected cross section is over an order of magnitude below this level 
of sensitivity. Stable sleptons can also be produced from cascade decays of 
heavier particles such as charginos. Such decays would also produce charged 
leptons and neutrinos, and the electron and I^t trigger data samples are 
searched to add sensitivity to these decays. The efficiency for this is very 
model dependent, and we quantify it only for a single point in the GMSB 
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parameter space which makes the three charged sleptons nearly degenerate 
with masses ~ 105 GeV/c^, slightly above the existing limits |^. The mod- 
ihed kinematics from the decays increases the efficiency to 6.7% for the muon 
trigger data set. Including the.^T and electron triggers increases it to 8.2%. 
The.^T trigger and isolation requirement introduce additional systematic un¬ 
certainties from the modelling of initial and hnal state radiation, making the 
total systematic uncertainty 12.5%. When cascade decays from all produc¬ 
tion modes are included, the cross-section limit is lowered to 550fb compared 
to the model prediction of 80fb. 

In summary, we have searched for long-lived charged massive particles 
in 90 pb“^ of data at CDF. No excess over background was observed. We 
derive cross-section limits using reference models for the two cases of strongly 
and weakly produced particles. In the strongly interacting case, these limits 
extend the excluded mass region to about 200 GeV/c^. 
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tions for SVX’ (open points) and CTC (solid points). The CTC prediction 
is a fit including detector effects. The SVX’ prediction is the Bethe-Bloch 
formula. The agreement is good in the low and high /d 7 regions important 
to this search. 
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Figure 2: Observed y^dE/dx distribution for tracks passing all the cuts for 
the strong production search in the muon trigger and If^T trigger data sam¬ 
ples [0. The curves are the expected background distributions which have 


an uncertainty of about 15% which is shown by the gray bands. 
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Figure 3; Limits set at 95% confidence level on the production cross sec¬ 
tion of long-lived fourth generation quarks are compared to the theoretical 
prediction. 
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